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We derive closed expressions for the nuclear breakup cross sections in the adiabatic limit using the 
Austern-Blair theory. These expressions are appropriate for the breakup of weakly bound nuclei. 
The concept of an exit doorway that mediates the coupling between the entrance channel and the 
breakup continuum is used. We prove the validity of the scaling law that dictates that the nuclear 
breakup cross section scales linearly with the radius of the target. We also compare our results for 
the nuclear breakup cross section of ^^Be, *B on several targets with recent CDCC calculation. 
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The breakup of nuclei is a common occurance when the 
bombarding energy is high enough and/or the binding en- 
ergies are sufficiently low. In the case of weakly bound 
nuclei the threshold for breakup is small and more so 
for bound unstable nuclei. The mechanism of breakup is 
assumed to consist of elongating the projectile, through 
the action of the interaction, which eventually leads to 
the production of two or more fragments. This interac- 
tion is composed of a short range, nuclear piece and a 
longer ranged electromagnetic one. A debate has been 
going on in the literature concerning the way the nuclear 
part of the breakup cross section depends on the mass 
of the target nucleus which supplies the interaction. In 
most references [J 0, it is assumed that the depen- 
dence goes as the cubic root of the mass number. In 
reference however, it is claimed that this dependence 
is more like linear! In a recent paper through a care- 
ful Continuum Discretized Coupled Channels (CDCC) 
calculation, the former dependence ( A^^^ ) has been es- 
tablished, which corroborates the contention that the nu- 
clear breakup cross section should follow the prediction 
of the Serber model . 

It is interesting to compare the numerical CDCC calcu- 
lation alluded to above with those of simpler analytical 
models. Specifically, the Ausern-Blair adiabatic theory 
for inelastic scattering comes to mind. If one assumes 
that the breakup proceeds through a so-called exit door- 
way 0, B , then the process can be treated as an in- 
elastic excitation. The idea of exit doorway has been 
used in the case of the influence of breakup on fusion 
0, [1] and in the excitation of giant resonances with 
success. In a very recent paper [10| a comparison of a pre- 
liminary CDCC calculation for the breakup cross section 
of the system ^He-|-^^A1 at low bombarding energies with 
a simple formula derived by us] using the Austern-Blair 
model, showed that such an idea is quite reasonable and 
encouraged us to pursue the matter further. We do this 
in the present paper, where we fully develop the Austern- 
Blair model for the nuclear elastic breakup reaction cross 
section assumed to proceed through the excitation of an 
exit doorway 0, H, 01 • 

The exit doorway concept has been used in the devol- 
opment of reaction theories invloving the the excitation 



of a doorway in the final state, in contrast to the con- 
ventional cases where such resonances are populated in 
the entrance channel In the breakup reactions of 

halo nuclei one may envisage that the process proceeds 
through the breakup doorway dipole, quadrupole etc.) 
into the continuum. As such, the detailed description of 
the exclusive reaction, where the final channels are espec- 
ified, will necessarily contain the full information about 
the exit doorway ( its energy, width etc.). This is the 
case that was encountered in the theory of the excitation 
of multiple giant resonances and of the influence of 
the pygmy resonance on the fusion of halo nuclei 7] . In 
the current paper we will be content with the inclusive 
quantity of the integrated breakup cross section and the 
only reference to the exit doorway is made implicitly as a 
final state that has to be populated for breakup to occur. 

The full Hamiltonian which describes the colliding ions 
can be written as 

H = Ho+F (1) 

where Hq = Iiq + K + V = ha + H^^'' is diagonal in open 
channel space, ho is the intrinsic part that describes the 
structure of the projectile and the target nuclei, K is the 
kinetic energy operator and V is the optical potential 
which contains the complex nuclear plus the Coulomb 
parts. The operator F describes the coupling among the 
open channels. 

The intrinsic Hamiltonian Hq, which for simplicity is 
taken here to represent the excitable projectile nucleus 
with the target considered structureless, is now written 
as: 

ho = \(j)o > Eq < (l)o\ + \d > Ed < d\ + > E, < i\ 

i 

+ Y^[\d > A, <i\ + \i> A* < d\] + 

i 

ij 

The first three terms on the RHS above refer to the 
ground, exit doorway and discretized continuum states. 
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respectively. The fourth term couples the doorway to 
the discretized continuum states, and the last term rep- 
resents the continuum-continuum coupling. If we remove 
the doorway from the above we get the Continuum Disc- 
trctized Coupled Channels ( CDCC) intrinsic Hamilto- 



ho = 100 > £^0 < 0ol + XI I* > < * 

i 

J2[\(j)o > A,<i\ + \i> A* < (j>o\ 



E 



[\i > f7y < j\+cc] 



(3) 



The exit doorway modulated CDCC Hamiltonian, Eq. 
(2), is our subject of study here. A full development of 
this new CDCC will be left for a future work. Here we 
concentrate our effort on understanding the consequence 
of reaching the breakup continuum from the entrance 
channel only through the exit doorway \d >. For this 
purpose we ignore the last term in Eq. (2) and remind 
ourselves that, whereas |0o > and \i > are eigenstates of 
/iQ, \d > is not. 

The full doorway-modulated CDCC equations can be 
obtained as follows. The full Schrodinger equation of the 
colliding system is, 

[E ~ (Ha + F)]\^j >^ (4) 

which when projecting onto the different channels gives: 



(5) 



iE~E,~Hl°^)ijt =F,o^P+ (6) 
We now invoke the exit-doorway hypothesis, 



Fqi — Foda*di Fm = Fdoa* 



(7) 



The overlaps adi and a*^ and can be easily obtained 
from Eq. (2) (without the last term)0, Q, 

= (rV2^)/[(£;. - Edf + (ri/2)2] (s) 

where F-^ the exit-doorway spreading width describ- 
ing its average coupling to the continuum states of the 
projectile, is related to the factors through. 



F^ =2^|AJ2 



P 



(9) 



where |AiP is an average value and p is the average den- 
sity of discretized continuum sates in the vicinity of d. 
Clearly the need to the continuum-continuum coupling 
terms would be very important if exclusive cross sections 
are to be calculated, since through them ( and through 
the doorway) the elastic channels coupling to the breakup 



channel continuum can be fully acounted for. Including 
the C-C coupling term, would result in a more compli- 
cated expression for \adi\'^ than that of EqEl 

Equations (5) and (6) can be recast into the following, 
after setting Eq = and = 0, 



^ad^Fod^. 



(+) 



(10) 



iE-E,^H<^'^)4+^=a*,Fdo^l;^ (11) 

The breakup cross section, within the exit doorway 
model then becomes, 



Kq 



M|<^(-)|i.,„|4+)>|2 (13) 



Where the sum over i has been performed by appropri- 
ate contour integration over Ei . Note that the Q-valuc in 
"0^ ^ is complex owing to the non-zero width of the exit 
dooway whose energy is Ed — «rjj/2. A simple way to see 
how the complex Q-value arises is to eliminate ip^^^ in 
Eq.(5) in favor of V'o^'' by employing Eq.(6), which gives 



^ ^ ^^(0) , - F,Q'ip^„'^\ With this Eq.(5) becomes 



the exit doorway hypothesis, the polarization potential 
contribution, J^i ^oi „ „ , . F^ becomes 



:i^zo)4^^ = 0. With 



V^/2ti 1 



{E, - EdY + (Fi/2)2 E-E,- Hi 



r(o) 



-Fdo 



Fod 



1 



E~{Ed~ iri/2) - H^"^ + ie 



Fdo 



This suggests defining the exit-dorway scattering wave 
function by setting H^^^ — H^'' such that Eqs.(lO) and 
(11) become: 



{E - {Ed - ^ry2) - Hf'^)^'^;'^ = FdQ^^Q 



',(+) 



(14) 



(15) 



The "inelastic" cross-section is thus given by Eq.(13) 
above with the aforementioned proviso that the Q-value 
of the excited state is complex. The width of this Q- 
value is a measure of the continuum contribution to the 
coupling. 

At this point we comment on the inclusion of the 
continuum-continuum coupling, namely the last term in 
Eq.(3). In this situation the amplitudes adi are obtained 
by matrix diagonalization and, among other things, the 
resulting overlap probability lofdip, deviates from the 
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Breit-Wigner form of Eq.(8). A possible form which may 
incorporate some of the c-c effects is a Lorentzian: 



The above form results in an equation for 4'^^'^ with a 
modified form factor which depends on the position and 
width of the exit doorway V^o ~ /(^-d, ri)Vdo where 
f{Ed,T\j) is generally complex. Accordingly the cross- 
section would be: a = \f{Ed,T\j)\'^aDWBA- In the hm- 
iting case of Vd <C Ed, the factor f{Ed,T\j) is approx- 
imately given by (1 — ^^Y^'^ resulting a cross-section 
pi 

given by: a~ {1+2e2^'^dwba- In the case of coupling to 
the breakup continuum considered here, the other limit, 
3> Ed is more appropriate, as Ed is roughly given by 
the Q- value of the breakup (< IMeV) while measures 
the extent in continuum excitation the discretization is 
performed 10 MeV). The function f{Ed,T^J can be 
calculated in such a situation, but we leave this for a fu- 
ture investigation. The important point we are making 
here, is that a DWBA calculation with complex excita- 
tion energy in the final state, and with a form factor of 
the type f{Ed,T'^)Vdo, should be an adequate candidate 
to treat the elastic breakup process. 

In the following we take the exit doorway to be excited 
sates of different multipolarities and use the Austern- 
Blair sudden/adiabatic theory [13, We employ the 
Distorted Wave Born Approximation for ip^'^'^ and "* ■ 

The elastic breakup cross section and its dependence 
on the target masss can be analysed within the Dis- 
torted Wave Born Approximation (DWBA). If we treat 
the breakup as an inelastic multipole process, the ampli- 



tude T, 



LM 



(+) 



> would look like: 



h 

< IfL: -MM\liQ > Rf-^d^ikf, h) 

g».,^(fe,)+.^(fe.)y^^ _^^(^^0). (16) 

The unpolarized cross section of the dipole transition 
is then obtained from the expression 



II k 



M=L 



LM\ 



(17) 



M=-L 



The radial integrals i?^ i.{kf,ki) for pure nuclear ex- 
citation are given by. 



Rf^^i^{kf,h) = {An/kfh) I drfij{kf,r)FL{r)fu{h,r) 

(18) 



Where the form factors FL{r) are given by the following 
expressions for the monopole,L = 0, dipole, L = I and 
quadrupole excitations L = 2 [l^ . 



dr 



Fi{r) 



,N).i^.l\Rp^^dV{r) , ^R^d''V{r 
-«i (o)( 



Rp 



dr 
F2(r) 



(N)dV{r) 



dr 



(19) 
(20) 
(21) 



with Ai?p = Rn — Rp being the difference between the 
rms radii of the neutron and proton distributions of the 
projectile, and V{r) is the elastic scattering channel op- 
tical potential. The quantities Rn and Rp can be ex- 
tracted from the analysis of Refs.flGl. \\^. In Ref. a 
power expansion in Ai?p was employed in the analysis 
of a-inelastic scattering from neutron skin nuclei. 

In the adiabatic limit, ki = kf = k, and for large or- 
bital angular momenta. If = li = I, the radial integral 
can then be evaluated in closed form following the pro- 
cedure of Austern and Blair [l3, [13] ■ For the dipole and 
quadrupole cases we have 



i?«(fc) 



[N) 



(k) 



dl 



Rfsj^. 

h' dp ' 



(22) 
(23) 



where 5^^^ is the nuclear deformation length given by 



: P^l'^Rp with being the nuclear deformation 



parameter and i?p is the radius of the excited projectile. 

The above expression for the radial integrals can be 
associated with the nuclear elastic breakup radial inte- 
gral. Thus we can obtain analytical expression for the 
integrated nuclear breakup cross section by simply inte- 
grating the cross section formula, eq.(2). In performing 
this calculation the angular momentum coupling coeffi- 
cients are evaluated exactly and the sum over li can be 
performed by putting the Coulomb phase shifts both as 
functions oilf = I. The amplitude of eq. (13) is given 
now by 



oo 



LM 



= ^V2^(2^+l)l/2i?i (fc)e2-,(fe)i^ _^^(^?^o). (24) 

1=0 



with the condition that Tlm = if i + Af is odd. The 
integrated pure nuclear breakup cross section containing 
dipole and quadrupole contributions then becomes the 
following 



i=0 



,dS'W(fc), 



dl 



(25) 
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where terms proportional to the second derivative of 
sl^\k) have been dropped. 

A simple estimate of the above formula can be made by 
approximating the sum in I by an integral in A = / + 1/2: 



E(2^ 

1=0 



1)1 



dl 



2A| 



dX 



(26) 

Assuming a real nuclear S-matrix which depends on 
A through [1 + exp(A — A)/A]~^ then the derivative of 
S would peak around the grazing angular momentum 
A with a width given by A. The integral (|26p is then 
obtained as: / = for A/A >> 1. Using A = kR 

A = ka, with R = ro{Ap^'^ + A^'^) and a being the 
diffuseness of the optical potential we find the simple 
formula for a. 



3a 



(27) 



where c is a constant normalization factor which depends 
among other things on the exit doorway nature of the 
excited state exemplified by the factor f{Ed,T [d)- It is 
clear that a depends linearly on the radius of the tar- 
get and, more importantly on the square of the nuclear 
dipole and quadrupole deformation lengths. Thus, the 



,1/3 



-dependence is established. 



In the calculation to follow we use the cluster model to 
calculate de deformation lengths for the different multi- 
polarities [l^, [13, [2lj This model assumes that the pro- 
jectile is composed of two clusters, a core of mass and 
charge ac and and a "valence" particle with Of, and zjj. 
The separation energy is denoted by Q, the Q- value of the 
breakup. Calling the spectroscopic factor of finding the 
cluster configuration in the ground state of the projectile, 
S one obtains the following expression for the distribu- 
tion of B{E\) in the excitation energy [T^ [2]|. 



rfSGBA) _ .2 



A-l 



(A!)2(2A + 1)( — )^x 

Q'/^{E.,~Q)^+'/^ 

p2A+2 ^ 



A^ 



^ fe^ (28) 



where A^o is normalization factor which takes into ac- 
count the finite range, rQ of the c+b potential. The latter 
is assumed to be such as to give a Yukawa type wave func- 
tion at large distances, ipbc{r) — Nq^JK /{2t:)^ — with 

K = \J 2^,bcQI'h^ and — --^^=. It is easy to obtain 
B{E\) by simply integrating of Eq 
the expression: 



(y + 1)2^+2 



dy = 



2A-I-3, 



2A-I-1 
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FIG. 1: CDCC calculations for the nuclear breakup (dots) 
compared to the results of Eq. (|27p . See text for details. 



We get for the cluster-model deformation lenghts Sf 



and $2 the following: 



2ttAp 



,B{E1) 



3 ZpNp' 
3 ,AcZh - 



SZpNp IGtt 



■(- 



(^V = ( 



Af 



47r 



Q 



(30) 



2,ZpRp' 



4:71 A f 



'SZpRp 327r \Hb, 



V ,AlZk 



and employing ^j^^ ^^^^ 



(2A + l)!sm[(A 



n (^+^fc) 



2 

(29) 



where p{= b + c) refers to the projectile. 

For our three nuclei discussed here, we have 
iiBe=i°Be+r7., **B='^Be-|-p and '^Be='*He-h^He, which 
define their cluster character, with the corresponding 
breakup Q-values, 0.504 MeV, 0.137 MeV and 1.587 
MeV. The factor N^S could be related to the Asymp- 
totic Normalization Coefficient (ANC) of the bound state 
wave function and is taken as a parameter to be adjusted 
so as to account for the experimentally known BlEX). 

Simple estimate of can be obtained from [ij, [l5| 

r ~ " a"^" ■ ^^^^^ present 
the results of the deformation lengths obtained using the 
cluster model of formulas [5D1 and [5T] For ^^Be, we used 
B{El) = 1.05 ± 0.06 e^fm^ Q and we get {5[^^f = 
0.71 ± 0.04 fm2. Further, (4^^ )2 = 1.27 ± 0.25 fm^ from 
the same reference. 
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TABLE I: Deformation lengths for the '^Be, **B and "Be pro- 
jectiles. The deformation lengths for ^Be and *B have been 
calculated using formulas (|30|l and (|3H) using A^oS = 1. For 
the ^^Be the Si and 82 axe the values from Ref. [l[ 



Projetile 


(5r')^(fm^) 






c 


■'Be 


0.41 


1.54 


0.00576 


1.37 


«B 


1.62 


3.24 


0.0179 


0.307 


"Be 


0.84(2) 


1.27(25) 


0.0214 


2.52 



In figure 1 we compare our results, Eq. ([?7)) using 
a = 0.65 fm with the CDCC calculation of Ref. [5] at 
Eiab = 200 MeV.A. Clearly we underestimate the CDCC 
calculation. The reason resides in the neglect, in our 
model, of the higher-order channel coupling terms al- 
luded to above. We also show in figure 1 the comparison 
with the CDCC calculation for and ^Be using the for- 
mula ip?)) . Clearly the scaling law is better obeyed in the 
"normal" nucleus ^Be as has already been discussed in 
Q . The value of the normalization c is close to unity for 
the "normal" nuclei ^Be For the *B the normalization 



is very close the average value of the Asymptotic Nor- 
malization Coefficient ANC = 0.45 measured in ref. [2^1 . 
For the ^^Be a higher normalization is obtained probably 
due to higher order effects which are not accounted by 
our DWBA description. 

In conclusion, we derived an expression for the nu- 
clear breakup cross-section using the Austern-Blair the- 
ory. The obtained cross-section exhibits the scaling law 
and should serve to supply a simple mean for an esti- 
mate of the nuclear breakup contribution. The expres- 
sion found should be contrasted with the purely geomet- 
ric Serber-like expression a — 



Acknowledgements 

This work is supported in part by Fundagao de Am- 
paro a Pesquisa do Estado de Sao Paulo (FAPESP) and 
the Conselho Nacional de Desenvolvimento Cienti'fico e 
Tecnologico (CNPq). M.S.H is the Martin Gutzwiller 
Fellow, 2007/2008. 



[1] N. Fukuda et al. Phys. Rev. C70(2004) 054606. 

[2] R. Palit et al. Phys. Rev. C68(2003) 0344318. 

[3] T. Aumann, Eur. Phys. J. A26(2005) 441 

[4] M. A. Nagarajin, C. H. Dasso, S. M. Lenzi and A. Vitturi, 

Phys. Lett. B, 503(2001) 65. 
[5] M. S. Hussein, R. Lichtenthaler, F. M. Nunes and I. J. 

Thompson, Phys. Lett., B640(2006) 91. 
[6] R. Serber, Phys. Rev. 72(1947) 1008. 
[7] M. S. Hussein and A. F. R. de Toledo Piza, Phys. Rev. 

Lett. 72(1994) 2693. 
[8] M. S. Hussein, M.P.Pato and A. F. R. de Toledo Piza, 

Phys. Rev. C52(1995) 846. 
[9] L.F. Canto, A. Romanelh, M. S. Hussein and A. F. R. de 

Toledo Piza, Phys. Rev. Lett. 72(1994) 2147. 
[10] E.J. Benjamin et al., Phys. Lett. B 647(2007) 30 
[11] N. Auerbach and V. Zelevinsky, Nucl. Phys. A781(2007) 

67. 



[12] N. Austern and J.S. Blair, Ann. Phys. (NY) 33(1965) 15. 

[13] W.E. Frahn, Nucl. Phys. A272(1976) 413. 

[14] G.R. Satchler, Nucl. Phys. A195(1972) 1. 

[15] G.R. Satchler, Nucl. Phys. A472(1987) 215 

[16] G. Alkhazov et al. Phys. Rev. Lett. 78(1997) 2313. 

[17] I. Tanihata et al., Phys. Lett. B289(1988) 261. 

[18] A. Krasznahorkay et al., Phys. Rev. Lett. 66(1991) 1287. 

[19] C.A. Bertulani. G. Baur and M.S. Hussein, Nucl. Phys. 

A526(1991) 751 
[20] C.A. Bertulani and A. Sustich, Phys. Rev. C46(1992) 

2340 

[21] C.A. Bertulani, M.S. Hussein and G. Miinzenberg, 

Physics of Radioactive Beams (Nova Science, New York, 
2001) 

[22] L. Trach, F. Carstoiu, C.A. Gaghardi and R.E. Tribble 
Phys. Rev. Lett. 87(2001) 271102-1 



